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We employ a straightforward relation between mutually unbiased and Bell bases to extend the latter in terms 
of a direct construction for the former. We analyze in detail the properties of these new generalized Bell states, 
showing that they constitute an appropriate tool for testing entanglement in bipartite multiqudit systems. 
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I. INTRODUCTION 

Entanglement is probably the most intriguing feature of the 
quantum world, the hallmark of correlations that delimits the 
boundary between classical and quantum behavior. Although 
some amazing aspects of this phenomenon were already no- 
ticed by Schrodinger in the early stages of quantum theory Ql, 
it was not until quite recently that it attracted a considerable 
attention as a crucial resource for quantum information pro- 
cessing J2]. 

The simplest instance of entanglement is most clearly il- 
lustrated by the maximally entangled states between a pair of 
qubits (known as Bell states), whose properties can be found 
in many textbooks Hfl. Despite their simplicity, they are of 
utmost importance for the analysis of many experiments fl]. 

In consequence, as any sound concept, Bell states deserve 
an appropriate generalization. However, this is a touchy busi- 
ness, since thoughtful notions for a pair of qubits, may be- 
come fuzzy for more complex systems. There are two sen- 
sible ways to proceed: the first, is to investigate multipartite 
entanglement of qubits. While the standard Bell basis defines 
(for pure states) a natural unit of entanglement, it has recently 
become clear that for qubits shared by moreparties there is a 
rich phenomenology of entangled states Jl uU, HUGH El- 

The second possibility involves examining bipartite entan- 
glement between two multidimensional systems |[T2l IT3I [Til 
Again there is no unique way of looking at the prob- 
lem, and different definitions focus on different aspects and 
capture different features of this quantum phenomenon. 

We wish to approach this subject from a new perspec- 
tive: our starting point is the notion of mutually unbi- 
ased bases (MUBs), which emerged in the seminal work of 
Schwinger 01711 and it has turned into a cornerstone of quan- 
tum information, mainly due to the elegant work of Wootters 

and coworkers SfHElUllHi. Since MUBs contain com- 
plete single-system information and Bells bases about bipar- 
tite entanglement, one is led to look for a relation between 
them. 

In this paper we confirm such a relation for qudits l23ll and 
take advantage of the well-established MUB machinery (in 
prime power dimensions) to propose a straightforward gener- 
alization of Bell states for any dimension. The resulting bases 
are analyzed in detail, paying special attention to their symme- 
try properties. In view of the results, we conclude that these 
states constitute an ideal instrument to analyze bipartite mul- 



tiqudit systems. 

H. BIPARTITE QUDIT SYSTEMS 

A. Mutually unbiased bases for qudits 

We start by considering a qudit, which lives in a Hilbert 
space Ttd, whose dimension d is assumed for now to be a 
prime number. The different outcomes of a maximal test con- 
stitute an orthogonal basis of Tid- One can also look for other 
orthogonal bases that, in addition, are "as different as possi- 
ble". 

To formalize this idea, we suppose we have a number of 
orthonormal bases described by vectors where £ (£ = 

0,1, . . . ,d — 1) labels the vectors in the nth basis. These are 
MUBs if each state of one basis gives rise to the same proba- 
bilities when measured with respect to other basis: 



1 



n 7^ n . 



Equivalently, this can be concisely reformulated as 



(2.1) 



(2.2) 



Note in passing that the Hermitian product of two MUBs is 
then a generalized Hadamard matrix, i.e., a unitary matrix 
whose entries all have the same absolute value 12411 . 

If one wants to determine the state of a system, given only 
a limited supply of copies, the optimal strategy is to perform 
measurements with respect to MUBs. They have also been 
used in cryptographic protocols l25ll . due to the complete un- 
certainty about the outcome of a measurement in some ba- 
sis after the preparation of the system in another, if the bases 
are mutually unbiased. MUBs are also important for quan- 
tum error correction codes ll26l, l27ll and in quantum game the- 
ory m Eh- 

The maximum number of MUBs can be at most d+ 1 l32ll . 
Actually, it is known that if d is prime or power of prime 
(which is precisely our case), the maximal number of MUBs 
can be achieved. 

Unbiasedness also applies to measurements: two nonde- 
generate tests are mutually unbiased if the bases formed by 
their eigenstates are MUBs. For example, the measurements 
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of the components of a qubit along x, y, and z axes are all un- 
biased. It is also obvious that for these finite quantum systems 
unbiasedness is tantamount of complementarity i33ll34tl . 

The construction of MUBs is closely related to the possibil- 
ity of finding of d + 1 disjoint classes, each one having d — 1 
commuting operators, so that the corresponding eigenstates 
form sets of MUBs Q5ll . Different explicit methods in prime 
power dimensions have been suggested in a number of recent 
papers JM, H3, [MS EH, but we follow here the one 
introduced in Ref. 042I1 . since it is especially germane for our 
purposes. 

First, we choose a computational basis \£) in Hd and intro- 
duce the basic operators 



X\t) = \t+1), 



Z\£)=u(l)\£) 



(23) 



where addition and multiplication must be understood modulo 
d and, for simplicity, we employ the notation 



=cu e = cxp(i27r£/d) 



(2.4) 



co = cxp(i27r / d) being a dth root of the unity. These opera- 
tors X and Z, which are generalizations of the Pauli matrices, 
were studied long ago by Weil [43]. They generate a group un- 
der multiplication known as the generalized Pauli group and 
obey ZX = to XZ, which is the finite-dimensional version of 
the Weyl form of the commutation relations [44]. 
We consider the following sets of operators: 



A(m)=X m , A(m,n) = Z m X r 



(2.5) 



with to = 1, . . . , d — 1 and n = 0, . . . , d — 1. They fulfill the 
pairwise orthogonality relations 



Tr[A(m) A^(to')] = d8 mm i . 
Tr[A(m, n) A 1 * (to', n')] = db mm * 5 nn > , 



(2.6) 



which indicate that, for every value of n, we generate a maxi- 
mal set of d — 1 commuting operators and that all these classes 
are disjoint. In addition, the common eigenstates of each class 
n form different sets of MUBs. 

If one recalls that the finite Fourier transform F is 14511 



a— i 
^ d t,t>=0 



then one easily verifies that 

Z = FX F^ 



(2.8) 



much in the spirit of the standard way of looking at comple- 
mentary variables in the infinite-dimensional Hilbert space: 
the position and momentum eigenstates are Fourier transform 
one of the other. 

The operators A(m, n) can be written as 



A(m, n) = e^C"'™) V n Z m U tn 



(2.9) 



where V turns out to be (d > 2) 

d-l 



1=0 



(2.10) 



and the phase 4>(m, n) is B46LI4" 

(f>(m,n) = uj(2~ 1 nm 2 ) . (2.11) 

Here 2 _1 denotes the multiplicative inverse of 2 modulo d 
[that is, 2" 1 = (<i+l)/2] and \£) is the conjugate basis, which 
is defined by the action of the Fourier transform on the com- 
putational basis, namely \l) = F\£). 

The case of qubits (d = 2) requires minor modifications: V 
is now 



V = 



1 ( 1 + i 1-i 
1-i 1 + i 



(2.12) 



while its action reads asVZV^ = —iZX. 

The operator V has quite an important property: its pow- 
ers generate MUBs when acting on the computational basis: 
indeed, if 



m = v n \i) 



(2.13) 



one can check by a direct calculation that the states \ip") 
fulfill (12.2b . which confirms the unbiasedness. If we denote 
Aw(m,n) = (£\A(m, n)\£'), according to Eq. j2.9i , we have 



A«,(TO,n)=e^ m <") mZ m \^) 



(2.14) 



Therefore, up to an unessential phase factor, A«/(to,ti) are 
the matrix elements of the powers of the diagonal operator 
Z in the corresponding MUB. This provides an elegant inter- 
pretation of these objects, which will play an essential role in 
what follows. 



B. Qudit Bell states 

For the case of two qudits^ a sensible generalization of Bell 
states was devised in Ref. j4g) , namely 

d-l 

|* mn ) = -= J2u(m£)\£) A \£ + n) B , (2.15) 
Vd l=0 

where, to simplify as much as possible the notation, we drop 
the subscript AB from \^ mn ), since we deal only with bipar- 
tite states. For further use, we also define 



"=0 



In the same vein, some gene ralized gates have been proposed 
to create these d 2 states l49l [50h . 
This set of states is orthonormal 

(2.17) 
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and allows for a resolution of the identity 

d-1 d-1 d-1 



m— 1 n— 



so they constitute a bona fide basis for any bipartite qudit sys- 
tem. As anticipated in the Introduction, there must be then a 
connection with MUBs. And this is indeed the case: it suffices 
to observe that the states (12.15b and ( 12.16b can be recast as 



I*, 



a— i 

= -nY. ^e(m,n)\e) A \e') B 



d-1 

-= ]T A u ,(m)\l) A \l') B 



(2.19) 



which can be checked by a direct calculation and A«/(m, n) 
and ku> (to) are the matrix elements of the operators ( 12.51 ). 

The matrices A possess quite an interesting symmetry prop- 
erty 

A u >{m, n) = Lu(m 2 n) A m (to, n) , A«<(m) = A^(to) . 

(2.20) 

In consequence, A(to) are always totally symmetric under 
the permutation of subsystems A and B and so are the cor- 
responding Bell states. Whenever oj(m 2 n) = ±1, A(m, n) 
are either symmetric or antisymmetric. This happens for 
mn = (mod d), and this is only possible for qubits: the 
symmetric matrices are A(0), A(l), and A(l, 0), while the an- 
tisymmetric is A(l, 1). The corresponding symmetric states 
are |* ) = 1*+). = 1*+). and |*i,o) = |$-), and 

l^i. l) = |^-) is the antisymmetric one. 

Finally, we can sum up the projectors of the bipartite states 
(12.15b over to, obtaining the following interesting novel prop- 
erty: 



d-1 



d-1 



m=0 



£ \*mn)(* m n\ = ® {X nl Z l ) B , 

e=o 

(2.21) 

d-1 1 d-1 

Y,\*m)(*m\ = 2 J2(X%®(X i ) B . 



rn=0 



i=0 



In words, this means that the sum of projectors over the index 
to is the sum of direct product of commuting operators for 
each particle. The proof of this statement involves a tedious 
yet direct calculation. 

For the case of two qubits, this implies that 

l*mlX*™l| = \[t + (XZ) A ® [XZ) B ] , 



m=0,l 



E \*™){* m \ = ^[l + (X) A ®{X) B ]. 



(2.22) 



m. BIPARTITE MULTIQUDIT SYSTEMS 
A. Mutually unbiased bases for n qudits 

The previous ideas can be extended for a system of n qu- 
dits. Instead of natural numbers, it is then convenient to use 
elements of the finite field F^n to label states, since then we 
can almost directly translate all the properties studied before 
for a single qudit. In the Appendix we briefly summarize the 
basic notions of finite fields needed to proceed. 

We denote as | A) (from here on, Greek letters will represent 
elements in the field F^™ ) an orthonormal basis in the Hilbert 
space of the quantum system. Operationally, the elements of 
the basis can be labelled by powers of the primitive element, 
which can be found as roots of a minimal irreducible polyno- 
mial of degree n over Z<j. 

The generators of the generalized Pauli group are now 



X fl \\) = \\ + I i), ^|A)= X (A/i)|A), 



(3.1) 



where xW i s an additive character (defined in the Ap- 
pendix). The Weyl form of the commutation relations reads 
as Z^X V = x{^v)X v Z^. 

In agreement with ( 12.5b , we introduce the set of monomials 



A(M) = x i 



ii ■■ 



A(/i, v) = Z„X. 



VfJ, 1 



(3.2) 



and their corresponding eigenstates also form a complete set 
ofcT + 1 MUBs. 

The finite Fourier transform now is [51] 



and thus 



f=4=E X(AA')|A)(A'|, (3.3) 



A, A' EE 



Zp = F X^ F' 



(3.4) 



The rotation operator V v transforms the diagonal Z^ into 
an arbitrary monomial according to 



A( ( ,,p) = e i ^^V v Z a VJ, 



rt 



(3.5) 



and is diagonal in the conjugate basis (defined, as before, via 
the Fourier transform |A) = F |A)) 



A 



ca,|A)(A| 



(3.6) 



where the coefficients c\ v satisfy the following relation 

CO;/ = 1 , Ca+A' v C*Xu = CAVX(-^A'A), (3.7) 

When d ^ 2, a particular solution of Eq. ( 13.7b is 

c A „ =x(-2" 1 A 2 ^). (3.8) 
Again, if we define the states 

W> = W>, (3-9) 

they are unbiased and Aaa' v ) are the matrix elements of 
the diagonal operator Z^ on the corresponding MUB 



m=0,l 



K X y{^ V ) = e^\rx\Z li \Vy) 



(3.10) 
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B. Multiqudit Bell states 

For a bipartite system of n qudits, it seems natural to extend 
the previous construction (12. 1 9b by introducing the d 2n states 

I*/-) = -i=E A ^'(/^m)^%, 



l=£A AA ,(/i)|A) A |A% 



(3.11) 



A, A' 



Accordingly, the associated Bell states are (apart from an 
unessential global phase) 



\ I fill * — * 



I*, 



* £|AU|A + ,) B 



(3.12) 



where tr, in small case, denotes the trace map in the field. 
Clearly, all the states l^o) an d l^u) are symmetric. The anti- 
symmetric states [i.e., A\y(fi, v) = —A\>\(fi, v)\ are defined 
by the pairs (jtt, v) such that 



tr(V/ii 2 ) = 1 



(3.17) 

Finally, a property similar to ( 12. 2U is fulfilled: summing 
up the projectors over /x one obtains 

E l*Mv)<*^l = J2^ X ^ Z -^ A ® (^A^a)b , 



(3.18) 



M A 

whose interpretation is otherwise the same as for qudits. 
C. Examples 



which look as quite a reasonable generalization. One can 
prove the orthogonality 



(3.13) 



and the completeness relation 



E 1*^X^1 +EI^)(*mI=1. (3-14) 

which confirms that they constitute a basis. Moreover, the 
reduced density matrices for both subsystems are completely 
random 



(and other equivalent equation with A and B interchanged) 
showing that they are maximally entangled states. 

The concept of symmetric and antisymmetric states can be 
worked out for systems of n qubits, which constitutes a non- 
trivial generalization of our previous discussion 1161 El. The 
symmetric states [i.e., Aaa'(Mi^) = A-A'aCMi correspond 
to those pairs (//, v) such that 



trO/r 2 ) = 0. 



Since we are dealing with n-qudit systems, we can map the 
abstract Hilbert space Ti^ into n single-qudit Hilbert spaces. 
This is achieved by expanding any field element in a conve- 
nient basis {8j} (with j = 1, . . . , n), so that 



A 



(3.19) 



where ^ G Z^. Then, we can represent the states as |A) = 
\£i, ... , and the coefficients ^ play the role of quantum 
numbers for each qudit. 

For example, for two qubits, the abstract state (|0) + 
|er 3 ))/\/2, where a is a primitive elements, can be mapped 
onto the physical state |00) + |10))/\/2 in the polynomial ba- 
sis {1, cr}, whereas in the selfdual basis {er, a 2 } it is associated 
with (| 00) + |11})/V2. Observe that, while the first state is 
factorizable, the other one is entangled. 

The use of the selfdual basis (or the almost selfdual, if the 
latter does not exist) is especially advantageous, since only 
then the Fourier transform and the basic operators factorize in 
terms of single-qudit analogues: 



(3.20) 



(3.16) 



For a bipartite 4x4 system the states are represented as 
| A) = \£i,£ 2 ) with £j E Z 2 . The Bell basis can be expressed 

as 



\mi, m; m 2 , n 2 ) = ^ '— E (-l) mifl+m2 ' 2 \h + mm 2 + m 2 n 1 J 2 + mmi + m 2 n 2 ) A \£i, £ 2 )b , 

(3.21) 
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The conditions 



mi(i 2 + m 2 ni 



0, 
1, 



(3.22) 



determine the symmetric and antisymmetric states, respec- 
tively. The solutions of this equation show that there are 10 
symmetric states and 6 antisymmetric ones, whose explicit 
form can be computed from previous formulas. 

Before ending, we wish to stress that so far we have been 
dealing with systems made of n qudits. However, sometimes 
they can be treated instead as a single 'particle' with d n levels. 
For example, a four-dimensional system can be taken as two 
qubits or as a ququart. If, for some physical reason, we choose 
for the quqart, we can still use Eq. ( 12.151 1, as in Ref. l48ll . even 
if now the dimension is not a prime number. However, if we 
proceed in this way the resulting basis contains 6 symmetric 
and 2 antisymmetric states, while the other 8 do not have any 
explicit symmetry, contrary to our results. 



IV. CONCLUDING REMARKS 

In summary, we have provided a complete MUB-based 
construction of Bell states that fulfills all the requirements 
needed for a good description of maximally entangled states 
of bipartite multiqudit systems. 

Mutually unbiasedness is a very deep concept arising from 
the exact formulation of complementarity. The deep connec- 
tion shown in this paper with Bell bases is more than a mere 
academic curiosity, for it is immediately applicable to a vari- 
ety of experiments involving qudit systems. 



APPENDIX A: FINITE FIELDS 

In this appendix we briefly recall the minimum background 
needed in this paper. The reader interested in more mathe- 
matical details is referred, e.g., to the excellent monograph by 
Lidl and Niederreiter lIBUl . 

A commutative ring is a nonempty set R furnished with two 
binary operations, called addition and multiplication, such 
that it is an Abelian group with respect the addition, and the 
multiplication is associative. Perhaps, the motivating example 
is the ring of integers Z, with the standard sum and multipli- 
cation. On the other hand, the simplest example of a finite 
ring is the set 1 n of integers modulo n, which has exactly n 
elements. 

A field F is a commutative ring with division, that is, such 
that does not equal 1 and all elements of F except have 
a multiplicative inverse (note that and 1 here stand for the 
identity elements for the addition and multiplication, respec- 
tively, which may differ from the familiar real numbers and 



1). Elements of a field form Abelian groups with respect to ad- 
dition and multiplication (in this latter case, the zero element 
is excluded). 

The characteristic of a finite field is the smallest integer d 
such that 



d 1 = 1 + 1 



1 = 



(Al) 



d times 



and it is always a prime number. Any finite field contains a 
prime subfield and has d n elements, where n is a natural 
number. Moreover, the finite field containing d n elements is 
unique and is called the Galois field F^™ . 

Let us denote as Zd [x] the ring of polynomials with coeffi- 
cients in Zd. Let P(x) be an irreducible polynomial of degree 
n (i.e., one that cannot be factorized over Z^). Then, the quo- 
tient space 7*d[X}/ P(x) provides an adequate representation 
of Fd" . Its elements can be written as polynomials that are 
defined modulo the irreducible polynomial P(x). The multi- 
plicative group of Fd» is cyclic and its generator is called a 
primitive element of the field. 

As a simple example of a nonprime field, we consider the 



polynomial x 



1 = 0, which is irreducible in Z2. If er is 



a root of this polynomial, the elements {0,1, ct, a 2 = a + 1 
a^ 1 } form the finite field F 2 2 and a is a primitive element. 
A basic map is the trace 



tr(A) = A + A 2 + . . . + A d " . 
It is always in the prime field and satisfies 
tr(A + A') = tr(A) + tr(A') . 
In terms of it we define the additive characters as 

~2ni 



(A2) 



(A3) 



(A4) 



X(A) = exp tr(A) 

L P 

which posses two important properties: 

X(A + A / )=X(A) X (A'), £ x(AA') = d%,A- 

(A5) 

Any finite field F^ can be also considered as an n- 
dimensional linear vector space. Given a basis {6j}, (j = 
1, . . , , n) in this vector space, any field element can be repre- 
sented as 



(A6) 



with tj £ Z,j. In this way, we map each element of F^™ onto 
an ordered set of natural numbers A (£1, . . . , £ n ). 

Two bases {81, . . . , 8 n } and {9[, . . . , 9' n } are dual when 



tr(M0=*M- 



(A7) 
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A basis that is dual to itself is called selfdual. 

There are several natural bases in F^™ . One is the polyno- 
mial basis, defined as 

{l^cr 2 ,...,^- 1 }, (A8) 

where a is a primitive element. An alternative is the normal 
basis, constituted of 

{a,^,...,^"'"}. (A9) 

The choice of the appropriate basis depends on the specific 
problem at hand. For example, in F22 the elements {er, <r 2 } 
are both roots of the irreducible polynomial. The polynomial 



basis is {1, a} and its dual is {a 2 , 1}, while the normal basis 
{a, a 2 } is selfdual. 

The selfdual basis exists if and only if either d is even or 
both n and d are odd. However for every prime power d n , 
there exists an almost selfdual basis of F<jn , which satisfies 
the properties: tr(0j0j) = when i ^ j and tr(#f ) = 1, with 
one possible exception. For instance, in the case of two qutrits 
F32, a selfdual basis does not exist and two elements {a 2 , cr 4 }, 
a being a root of the irreducible polynomial x 2 + x + 2 = 0, 
form a self dual basis 

tr(cr V) = 1 , tr(cr 4 cr 4 ) = 2 , tr(<7 2 cr 4 ) = . (A10) 
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